Testing gravity at the Second post-Newtonian level through gravitational deflection of 

massive particles 
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Expression for second post-Newtonian level gravitational deflection angle of massive particles is ob- 
tained in a model independent framework. Ambiguities associated with massless particle deflection 
can be avoided by using massive particles. Comparison of theoretical values with the observationally 
constructed values of post-Newtonian parameters for massive particles offers the future possibility 
of testing at that level competing gravitational theories as well as the equivalence principle. 
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I. INTRODUCTION 

A consequence of general relativity (GR) is that light 
rays are deflected by gravity. Historically, observations 
of this aspect of gravity provided the first proof in favour 
of GR. The effect is nowadays routinely used as a tool 
to study various features of the universe, such as view- 
ing fainter or distant sources, estimating the masses of 
galaxies etc. [1]. 

Like photons, particles having masses are also deflected 
by gravity. The general relativistic corrections in the 
equation of motion of a massive test particle moving in 
bound orbits have been studied with great accuracy in 
the literature [2]. These studies have great relevance 
in comparing theory with observations of gravitational 
waves from compact binary systems. However, for un- 
bound orbits of massive particle, the question of accuracy 
beyond first order has not received sufficient attention 
as of now. The main reason could be that the obser- 
vational aspect of unbounded massive particles was not 
very practical: there was no known astrophysical source 
of free point particles that can be detected easily with 
good angular precision. The situation seems to have im- 
proved somewhat. Recent theoretical studies [3] favor 
the existence of local astrophysical sources of relativistic 
neutral particles like neutrons and neutrinos with observ- 
able fluxes. Besides, high energy neutrons are produced 
during solar flares [4]. Moreover, with the advent of new 
technology new experiments have been proposed [5] , pri- 
marily to study gravitational deflection of light with high 
precision, in which laser interferometry will be employed 
between two micro-spacecrafts whose line of sight pass 
close to the sun. Hence there might be a possibility that 
in the future, neutron or some other neutral particle may 
be used in a similar experiment instead of photon thus 
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providing an opportunity of studying gravitational de- 
flection of massive particles. Henceforth, we use the ab- 
breviation for post-Newtonian as PN such that first-PN 
effect is of the order of (1/p), second-PN effect is of the 
order of (l//0 2 ) and so on. 



The expected angular precision of the planned astro- 
metric missions using optical interferometry is at the level 
of microarcseconds (/i arcsec) and hence these experi- 
ments would measure the effects of gravity on light at 
the second-PN order (c -4 ). Though measurements with 
massive particles at the level of microarcsecond accuracy 
is way beyond the present technical capability, it can still 
be cautiously hoped that astrometric missions in the dis- 
tant future using massive particle interferometry would 
have angular precisions close to that to be obtained using 
laser (optical) interferometry. Whatever be the technical 
scenario, a study of theoretical aspects of gravitational 
deflection of massive particles at the second-PN approxi- 
mation is useful in its own right. (To our knowledge, the 
deflection angle for massive particles has been theoreti- 
cally estimated in the literature with an accuracy of only 
first-PN order [6] so far). 



Accordingly, in the present article, the second-PN con- 
tribution to the gravitational deflection of massive parti- 
cles by a gravitating object will be estimated in a model 
independent way but with a special emphasis on the sun 
as gravitating object. The corresponding PN parameters 
for light deflection then follow as a corollary. The im- 
portance of this investigation is manifold: It allows us to 
obviate the ambiguities related with photon deflection in 
that order; it allows us to construct the true coordinate 
solar radius from measurements and moreover, it consti- 
tutes a possible futher test of the equivalence principle. 
These are discussed at the end. 
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II. GRAVITATIONAL DEFLECTION OF 
MASSIVE PARTICLES AT SECOND 
POST-NEWTONIAN ORDER 

We consider the general static and spherically sym- 
metric spacetime in isotropic coordinate which is given 
by (we use the geometrized units i.e. G = 1, c = 1) 

ds 2 = -B(p)dt 2 + A(p) (dp 2 + p 2 d0 2 + p 2 sin 2 0d<j) 2 ) (1) 

Restricting to orbits in the equatorial plane (9 = tt/2), 
the expression for the deflection angle for particles mov- 
ing with a velocity V as measured by an asymptotic rest 
observer can be written as [7] 



a (Po) = I{Po) - 7T 



(2) 



with (see appendix) 

I(Po 
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where p Q being the distance of the closest approach, 



(3) 



J = Po [A{p ){B-\p )-E)] 



1/2 



and 



E = l-V 2 



(4) 



(5) 



The PN formalism in some orders [7,8] is usually cm- 
ployed to describe the gravitational theories in the solar 
system and also to compare predictions of GR with the 
results predicted by an alternative metric theory of grav- 
ity. This method actually is an approximation for obtain- 
ing the dynamics of a particle (in a weak gravitational 
field under the influence of a slowly moving gravitational 
source) to one higher order in — (M is the mass of the 
static gravitating object) than given by the Newtonian 
mechanics. Following the PN expansion method, we as- 
sume the metric tensor is equal to the Minkowski tensor 
f\\u> plus corrections in the form of expansions in pow- 
ers of and considering up to the second-PN correction 
terms, we have 



M M 2 
B(p) = 1-2— +2A — 
P P 2 
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(6) 



(7) 



/3»,7» are the PN parameters (also known as the Ed- 
dington parameters), Si and can be considered as the 
sccond-PN parameters, i stands for either 7 or m denot- 
ing photons or massive particles respectively. Several of 
these parameters are different for different theories. In 
GR, all of them are equal to 1 as can be readily checked 
by expanding the Schwarzschild metric. 



We should note that the metric coefficients above are 
independent of any specific model; they result solely from 
the assumption of central symmetry. Starting with the 
expansion (6) and (7) per se, the expression for the angle 
of deflection for unbound particles up to the second-PN 
order follows from Eq.(3) and when — << V 2 , it works 
out to 



where 
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(8) 



(9) 



b m = ^ + (2 + 2 7m - M ~2 ( 7 m + ^r) (10) 

The above expressions are also valid for massless particles 
(a 7 , a 7 , & 7 ) as may be seen under the substitution V = 1. 
Clearly the deflection angle would be larger for particles 
in comparison to that of photons. Our calculation shows 
that the term representing the second order effect (b m ) 
contains only the three parameters /3 m , j m , S m and does 
not contain e m , a cubic order contribution. This implies 
that calculation of the deflection of unbound particle or- 
bits (including photons) by gravity to any given order 
needs only the knowledge of every term to that order in 
the expansions. In other words, to second-PN order, one 
needs to consider both in g oa and gij terms only up to 
Similarly, to third-PN order, which is not our inter- 
est here, we would need expansions of both the metric 
components up to order -^3- and so on. This is in con- 
trast to the case of planetary dynamics (bound orbits) 
where the calculation typically requires knowledge of g oa 
more accurately than (For instance, to calculate the 
planetary precession to the order of M , one expands g OQ 



for 



up to 2/3 ^ M while g^ is expanded up to only 27 ^ M ; 
next order accuracy, one would need to consider the com- 
plete expansion as given in Eqs.(6), (7) above so that the 
parameters S m , e m become important in this case). The 
deflection angle a m for the Schwarzschild spacetime can 
be obtained by taking (3 m — 7™ = S m = 1. 

The deflection angle also can be expressed in terms 
of coordinate independent variables, such as the impact 
parameter b which is the perpendicular distance from the 
center of the gravitating object to the tangent to the 
geodesic at the closest approach. In that case, p has 
to be replaced by b in Eq.(8), the Eq.(9) would remain 
unaltered but Eq.(10) would change to 
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V 2 



(11) 

Since impact parameter is the ratio of the angular mo- 
mentum and energy of the particle as measured by an 
observer at rest far from the gravitating object, it is a 
formally measurable quantity but is not very suitable for 
practical measurements [9]. 
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III. OTHER SIGNIFICANT EFFECTS 

In the present work, the mass distribution of the grav- 
itating object is assumed to be mainly spherically sym- 
metric; any deviation from such symmetry would produce 
their effects. The effect of quadrupole moment of the 
mass distribution on the deflection angle is proportional 
to J( 3 li f i R ; where R is the average radius of sun. Thus, 

Po 

even a small quadrupole moment parameter Jq could 
produce significant contribution to deflection. However, 
the effect is limited largely to the first-PN order (~ 0.1 
fi arcsec) while in the second-PN order the effect is too 
small (~ 1 0~ 7 — 10~ 8 p arcsec). If the gravitating object 
also has angular momentum, its effect on the deflection 
angle contributes to the second-PN order but it can be 
separated out. All these are discussed below. 



Effect of quadrupole moment of the mass 
distribution 



Theoretical value of solar quadrupole moment Jq, 
though it depends strongly on solar model used, is very 
small, of the order of 10~ 7 [5,8]. Since our study is aimed 
at sun as the gravitating object we have ignored higher 
order terms involving Jq. Thus due to the quadrupole 
moment of the mass distribution the effective mass pa- 



rameter becomes M, 



eff 



M 



1 



(3cos 2 6 



1 



which leads to the following corrections in the compo- 
nents of the metric tensors [5,8,10]: 



MR 2 

Sgoo(p) = Jq — 5— (3cos 2 9 - l) 



(12) 



order terms in 1/ p 2 containing Jq here because they have 
magnitude in the range 10~ 7 to 10~ 8 p arcsec, too small 
to be of any practical significance. We can justifiably 
ignore these second-PN quadrupole contribution. The 
quadrupole contribution to the deflection of light is given 
in Ref.[ll]. 



B. Effect of rotation 

The angular momentum of the gravitating object is 
assumed small as in the case of sun. The resulting leading 
term of the relevant metric tensor is 



AMa 



(15) 



where a is the angular momentum per unit mass of the 
object. The contribution of the rotation to the deflection 
angle is then given by [6] 



OLrot 



AMaV 
~1T 



(16) 



and 



The value of a can be positive or negative depending on 
the direction of rotation. When the angular momentum 
of the gravitating object is antiparallel with the direction 
of the incoming particle, a is positive and hence rotation 
causes larger deflection whereas for parallel angular mo- 
mentum, a is negative and the deflection angle will be 
less. Thus the rotational effect can be easily separated 
out from other contributions by studying the deflection of 
particles at two opposite sides of the gravitating object. 

The gravitational deflection angle of light with an accu- 
racy up to second-PN order readily follows from Eqs.(8)- 
(10), (14) and (16) using V = 1. 



MR 



s 9jk{p) = SjkliJQ — 5— (3cos 2 9 - l) (13) 



IV. DISCUSSION 



where R is the average radius of the mass distribution and 
9 is the angle between radius vector and the z-axis and 
hence in the equatorial plane 9 = tt/2. In the equatorial 
plane, the deflection caused by the quadrupole moment 
calculates to 



CiQM 



2J Q MR 2 



In 



V 2 J 



(14) 



Assuming R 
and taking V — 



p at the closest approach to the sun 
0.75, (3 m = -f m = S m ~ 1, for sun — 



2.12 x 10~ 6 , this first-PN quadrupole term cxqm ~ 0.1 p 
arcsec. It is roughly 7 orders of magnitude less than the 
first-PN deflection a m y- ~ 0.8 sec and is more than one 
order of magnitude less than the second-PN contribution 

b m (jf^j ~ 3.4 /i arcsec. We have not displayed the next 
higher order quadrupole terms involving Jq and second 



The study of gravitational deflection of massive parti- 
cles is important for several reasons which are discussed 
below. 

First of all, observations of gravitational deflection of 
massive particles with fj, arcsec precision could probe 
the gravitational theories at the second-PN level without 
any ambiguity. The second order predictions of gravita- 
tional deflection of light as evolved from different stud- 
ies are found to be ambiguous. This is because of er- 
roneous identification of theoretical deflection variables 
with the observables or measured quantities. For in- 
stance, when light ray just grazes the limb of the Sun 
and the isotropic radial coordinate p Q is identified with 
the measured (under Euclidean approximation) solar ra- 
dius, it gives a second-PN contribution of <~ 3.5 p arcsec 
to deflection angle in GR [12]. But if one uses the stan- 
dard Schwarzschild coordinates instead, the second-PN 
contribution to the deflection angle of light in GR fol- 
lowing from Eqs.(8)-(10) would be [^f - l] ^ which 
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is numerically about 7 p arcsec for light ray grazing the 
limb of the sun, provided the distance of closest approach 
in such coordinates is identified with the measured radius 
of the sun. The deflection angle can also be expressed in 
terms of coordinate independent variables, such as the 
impact parameter b. In that case, the second-PN contri- 
bution to deflection angle in GR becomes 7^"^- and 
when at closest approach, b is identified as the measured 
solar radius for light grazing the sun, the magnitude of 
second-PN deflection angle is ~ 11 \x arcsec [13]. Thus 
there exists a great confusion about the prediction of 
GR (or in fact of any viable gravitational theory) at the 
second-PN order. 

The fundamental reason for such an ambiguity is that 
the measurements of solar radius usually employ Eu- 
clidean geometry as an approximation [6] whereas the an- 
gle of gravitational deflection or other GR effects are prin- 
cipally based on the consideration of curved spacctime. 
But, comparing points in two different geometries i.e., in 
curved spacetime and flat spacetime is totally meaning- 
less [14]. The Euclidean approximation works tolerably 
well only up to the first order, that is, in weak field grav- 
ity caused by a source like Sun. The magnitude of the 
second order contribution is, however, of the same order 
as the error that arises due to such an approximation. 
Hence the numerical value of gravitational deflection an- 
gle of light can not be unambiguously predicted at the 
level of second-PN order within the theoretical scheme 
in vogue. Since the deflection angle for massive particles 
depends also on the velocity of the particle, the stated 
ambiguity can be easily avoided by measuring deflection 
angles for two or more velocities of the probing massive 
particle. 

In GR, a coordinate length like p a is not directly mea- 
surable, it can only can be indirectly "constructed" from 
the values of actual measurements. The PN parameters 
and also the otherwise unknown coordinate solar radius 
p (or equivalently, r D in standard Schwarzschild coor- 
dinates) can be constructed through least square fitting 
with the measured deflection angles a m and probing ve- 
locities V using the Eqs. (8)-(10). The idea is that the 
values of coordinates, p a and r Ql which refer to the same 
radial point, should be treated more like other PN pa- 
rameters (Pm, 7m, o~m) due to the fact that the "flat ge- 
ometry" spacetime points can not be algebraically identi- 
fied in a curved spactime [14]. Technically, however, the 
flat radial distances can be constructed by using met- 
ric gravity itself (Eddington expansion) in terms of a 
large set of unknown PN parameters (/3 m , 7 m , S m ) in- 
cluding p by fitting them with the observed data [7]. 
This method has been adopted, for example, by Shapiro 
and his group in the radar echo delay observations [15, 
16]. The resulting parameter values can then be com- 
pared with the theoretical predictions of deflection in GR 
as well as in other competing theories (like Brans-Dicke 
theory) in the second-PN order involving both massive 
and massless partcles. 

The study of gravitational deflection of massive parti- 



cles is also important in the context of testing the weak 
equivalence principle which is one of the fundamental 
postulates of general relativity. The principle states that 
the trajectory of a freely falling object is independent of 
its internal structure and composition. In other words 
all particles are coupled with spacetime geometry uni- 
versally. The principle has been tested with great accu- 
racy through different experiments, notable among them 
are the Eotvos type experiments [10] where comparison 
of gravitational and inertial masses of objects are made 
by measuring their acceleration in a known gravitational 
field. For massless particles like photons, however, such 
measurements obviously can not be performed. Instead, 
in such situation the principle is tested by examining 
whether the gravitational (second-PN) coupling param- 
eter 7 is universal for all particles, massive or massless. 
On the basis of supernova 1987 neutrino and optical data 
[17], a limit of | 7 7 — j m \< 3.4 x 10~ 3 has actually been 
found [18]. However, the mass of a neutrino m„ e is very 
small (if not zero); the present upper limit being m„ e < 3 
eV. Hence, a more conclusive experiment would be to ex- 
amine whether the gravitational couplings for photon and 
massive particles (other than neutrinos) are same or not. 
The observational value of | 7 7 — 7 r „ | should provide a 
direct answer as to the degree of validity of the principle 
in question. 

Remarks on particle deflection experiment 

The main concern, which is still far from resolved, 
is whether realistic experiments for observing gravita- 
tional deflection of massive particles can be devised or 
not. Here, we only speculate on some possibilities. The 
most important requisite in this context is to generate 
a beam of suitable test particles. Charged particles like 
protons or electrons have to be excluded as test particles 
because they suffer electromagnetic interactions by the 
interplanetary magnetic field. Among neutral particles, 
neutrinos are unlikely to serve the purpose as their speeds 
are almost, if not exactly, the same as the speed of light. 
Thus, neutrons seem to be the only feasible candidate. 
They are known to be produced during solar flares but 
they can at best be used to study the gravitational deflec- 
tion by an intermediate planet. If astrophysical sources of 
neutrons other than the sun are detected in future exper- 
iments, the problem of searching the test particle beam 
would be resolved automatically. Otherwise, one might 
hope to generate the beam only artificially. However, 
since neutrons are unstable with a mean lifetime of 886 
sec, only neutrons with a minimum speed of 0.75c can be 
used as test particles so that they do not decay during the 
travel from one micro-spacecraft to another. Though in 
(man-made) accelerator experiments (at earth) neutrons 
can be accelerated to such speeds, it seems improbable at 
the present stage of technology that neutrons can be ac- 
celerated to such high energies from a micro-spacecraft. 
This is a challenge for the future. 

Alternatively, stable massive objects, such as a bullet, 
can also be used as test particles but they must have 
a minimum speed of ~ 6 x 10 cm sec" 1 so that its to- 
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tal energy remains positive throughout the path (from 
one micro-spacecraft /earth to another spacecraft) and 
would not be captured by solar gravity. The fastest man- 
made object (Helios 2 solar probe) has a speed of about 
7 x 10 6 cmsec _1 . However, as already mentioned, accel- 
erating a material object to speeds of ~ 10 8 cmsec _1 and 
achieving the required level of p arcsec accuracy in the 
deflection angle is completely beyond present technical 
feasibility. 
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APPENDIX A 

The standard equations for a geodesic, namely 



A dx^ dx v 
~ds lr + 1 ds ds 

for the general metric (1) become 



= 



d 2 p A (dp 
ds 2+ 2A\ds 



~P 1 




= 



dH_ 




smttcosb 



(Al) 



(A2) 



(A3) 



d 2 ^ 
ds 2 



2 A \ dpdO „ d6d<p „ 
— I — t- ] -T--J — h 2cot9—-j- = 
p A I ds ds ds ds 



d 2 t B dp dt 
ds 2 B ds ds 



(A4) 



(A5) 



(primes denoting differentiation with respect to p) . If we 
choose 6 = 7r/2 and d9/ds = initially, Eq.(A3) warrants 
that they would remain the same always. Thus normal- 
izing time coordinate suitably, one obtains for orbits in 
the equatorial plane from Eq. ( A5) 



dt 
ds 



B- 1 



Integrating Eq.(A4) 



P 2 ^ = A-\J 2 
ds 



(A6) 



(A7) 



where J is a constant of integration. From Eqs.(A2), 
(A6) and (A7), one finally obtains 



Ap 2 J 2 \B 



17(37] + — --(o-£)=0 (AM 



which leads to Eq.(3). J can be conveniently expressed 
in terms of distance at closest approach. At the point of 
closest approach, dp/d<fi vanishes. Using this in Eq.(A8), 
one recovers Eq.(4). 
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